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I. INTRODUCTION 



In Ref. DU (hereafter referred as I) we have obtained a generalized and exact 
form of the master and Langevin equations. In this work we consider a partic- 
ular soluble model and numerical results related to such equations. We show 
that from the microscopic quantum-mechanical laws a sort of "irreversible" 
behavior emerges under the following conditions: a privileged initial condition 
(no correlation between the Brownian oscillator and the bath), a relevance 
criterion arising from the form of the Hamiltonian, which distinguishes the 
Brownian particle from the bath, and a natural extra-dynamical hypothesis in 
order to interpret the exact dynamical evolution. That is, to consider such a 
process during a time scale of observation smaller than the recurrence time and 
with a minimum resolution such that fluctuations cannot be seen (see figures 
of Sec. V). 

In Sec. II we exactly solve the particular model, while in Sees. Ill and 
IV we use this solution in order to obtain the coefficients of the master and 
Langevin equations. Sec. V contains the numerical results. 

II. THE EXACT SOLUTION OF THE MODEL 

Let h be the following Hamiltonian, 

N N 

h = n\Q) (n\ + £>„ K> M + E (\ n ) Kl + K> ■ W 

n=l n=l 

The eigenvalue problem, h \a u ) = a u \a u ) , is reduced to the algebraic system 

N 

Q (Q\a v ) + E 9n (unM = ot u (Q\a v ) , 

n=l 

(2) 

9n (£l\a v ) + u n (uj n \a u ) = a u (uj n \a u ) . 
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From the second of Eqs. (0) we have 

(uj n \a u ) = — — — (Q\a u ) , (3) 
for a v 7^ u n , W, n. Replacing Eq. (|3|) in the first equation of (||) we obtain 

N 2 

the secular equation a v — Vt = J2 —^z — • Finally, by using Eq. (||) and the 

n=l " n 
N 

completeness relation (Q\ + J2 I^Vi) (^nl = !■> we have 



n=l 



m\a„)\ 2 = ^ (1) 



l+ E 

n=l 

where we have taken into account that (a u \a u ) = 1. 

III. THE EXACT MASTER EQUATION 

As was showed in I the general solution of the master equation is given by 
(n = 0,...,N) (N n (t)) = J2 n Pnm(t) (N m (0)) , which in the particular model 
given by Eq. ([[]) reads 



N 



(N n (t)) = Pan(t) (N a (0)) + £ Pn n (t) (N n (0)) , 



n=l 



A' 



(5) 



(N n (t)) = P nQ (t) (N n (0)) + ]T Pnmit) (N m (0)) 



m=l 



From Eq. (13) of I we have 

N 

Pnm(t) = e~ i{a »~ av)t (a^} (a„\lp n ) {lJ) m M (*PmK) . 

(j,,v=0 

By taking into account the exact solution of the model [Eqs. (^) and (PS)] these 
probabilities are given by 

p fA ^ cos [(a„ - a v ) t] * 1 
Pnn{t) = 2. 2^ TTTTi — m-2 i /r >i — TF2 + 



^ KQla^r 2 \(Q\a v )\~ 2 ^ol(^l«,)| 
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N 9l 

+ 5> \{tt\a u }\~ 4 (a v - w n ) 2 ' ^ 
p ( t ) = 2 V 5 C ° S ^ ~ — ~ 

\x>v 

N a 2 a 2 

_|_ iJnilm 

where {|^ > , |^„)} = {|^> , K>} (n = 1, N). 

IV. THE EXACT LANGEVIN EQUATION 

As it was shown in I the solution of the Langevin equation is 
X{t) = \ [Ann(t) + A* Qn (t)} X(0) + 1 [A nn (t) - ^(t)] ^ + /(f), (7) 
which together with A nn (t) = £^ =0 \(n~\Z)r' 2 yields 

The coefficients result to be 

„ cos[(a 1 ,-a M )f]q^Q^ _ sin[Q„--a M )t]ag 

^ ' cos[(a^— Qjjjtjai/ ' ' cos[(a 1/ — a M )t]a 1/ ' v / 

These coefficients are independent of the initial condition which reflects the 
purely dynamical evolution of the Brownian particle. The privileged initial 
condition is only needed for computing the mean values of the relevant ob- 
servables leading to the apparently irreversible behavior from a macroscopical 
point of view. 

V. NUMERICAL RESULTS 



The parameters and functions that play the game are chosen as follows. 
The mean values are taken with respect to the canonical distribution for the 
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bath, which is in thermal equilibrium at temperature at t = 0. The initial 
population for the Brownian oscillator is fixed to the unity and its natural 
frequency Q = (3~ l . The frequencies of the bath oscillators are spaced with a 
constant step A and centered around Vt according to uo n = + A (n — , 
where A = uj n+ i —u n ;n=l,...,N. The coupling function is a Lorentzian given 
by o„ = , Z 4 " 2 ^'i , where a = A ( N ~ 2 ) , A more detailed discussion about this 
choice can be found in Ref. [§]. 

All the figures below correspond to N = 100, A = 0.018, = 1. 

A. Master equation 

In Fig. 1 we have the time evolution of (Nq) [given in Eq. (j5|)] departing 
from an initial population (Nq(0)) = 1. We can see that after a short non- 
exponential regime (known as Zeno period [Q]) the decay profile fits a decreas- 
ing exponential until it reaches the equilibrium value (Nn(t ^> fi _1 ,7 -1 )) = 
^ e /3n — 1^ =0.58 (where the Brownian oscillator thermalizes with the bath). 
It is the asymptotic value reached by (Nq) for times smaller than the recurrence 
time t r ~ . -, 2w r = 37,311, which is several orders of magnitude greater 

7 vam.(a v +i—a v ) ' ' to to 

than the oscillator period tq = 2tc/Q = 2tt (see a detailed explanation in Ref. 
@). Fig. 2 shows the contributions to (Nq) stemming from the survival prob- 
ability Pfifi(i) = (p* e~ tht Cl^ of the initially prepared unstable state \Q) , 
which reaches for t < t r a vanishing asymptotic value, and the contribution 
coming from the bath J2 n Pn n {t) (N n (0)) , which provides the equilibrium value 

to (Nq) . 

B. Langevin equation 



In Fig. 3 we show the damped oscillations of the Brownian particle com- 
paring it with the reconstruction and decay of (Nq) in the peak around the 
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recurrence time. Fig. 4 gives the time behavior of the "damping" coefficient 
T(t) re-scaled with respect to its asymptotic value 7 obtained in I through 
second order perturbation theory. It can be shown that T(t) essentially follows 
—2Rej^ at all times. We firstly see that T(t) grows from to 7, correspond- 
ing to the Zeno regime. Secondly, it oscillates around T = 7 as a consequence 
of the presence of fluctuations which modulates the exponential decay. After 
that, the amplitude of these oscillations increases in time because of the rel- 
ative variations between fluctuations and the exponential decay (— 2Re^^-) 
become more important when the system reaches equilibrium. Finally, we see 
the effect of recurrences also in T(t). It is interesting to stress the fact that 
F(t) allows us to visualize the presence of fluctuations which, at the same time 
scale, are hidden in the observable macroscopic profiles. 
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Figure captions: 

1. Population of the Brownian oscillator vs. t 

2. Survival probability vs. t 

3. Mean position of the Brownian oscillator vs. t 

4. Damping factor of the Langevin equation vs. t 
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Figure 1 
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Figure 3 
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Figure 4 
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